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I.  INTRODUCTION 


A.   BACKGROUND  INFORMATION 

As  the  size,  complexity,  and  hence  the  cost  of  modern 
warships  have  grown  in  the  last  thirty  years,  advances  in 
basic  research  have  been  made  and  the  computer  has  become  a 
useful  tool  to  the  researcher  and  designer.   The  application 
of  up-to-date  research  and  computer  technology  to  the  design 
of  modern  warships  is  common  but  there  remain  major,  costly 
components  relatively  unchanged  and  unreviewed  in  decades. 
It  is  in  the  interest  of  the  Navy  to  ferret  out  these  compo- 
nents and  re-examine  them. 

One  such  component  is  the  naval  condenser.   The  basic 
engineering  developments  currently  used  to  design  naval  con- 
densers are  from  the  Heat  Exchanger  Institute  (H.E.I.) 
Standards  for  Steam  Surface  Condensers  [1]  and  the  standards 
of  the  Tubular  Exchanger  Manufacturers  Association  (T.E.M.A.) 
[2].   The  Navy  has  further  utilized  the  H.E.I,  information 
in  the  Design  Data  Sheet  (D.D.S.)  [3].   While  these  standards 
have  proven  reliable,  their  performance  margins  and  extensive 
use  of  averaged  values  may  result  in  a  larger,  costlier  con- 
denser at  no  increase  in  performance  and  with  an  unnecessary 
margin  of  reliability.   These  performance  margins,  once 
necessary  to  insure  reliability,  are  now  called  into  question 
as  the  mysteries  of  condensation  are  unravelled. 


It  is  for  the  purpose  of  reducing  the  size  of  the  con- 
denser that  many  researchers  are  trying  to  increase  the  heat 
transfer  rate  per  unit  of  condensing  surface  area. 


B.   CONDENSATION  IN  A  NAVAL  CONDENSER 

A  typical  naval  condenser  uses  the  heat  transfer  mecha- 
nism of  filmwise  condensation  to  the  outside  of  horizontal, 
cylindrical  tubes.   Dropwise  condensation  is  therefore  not 
considered  in  this  thesis. 

The  inside  of  the  tube  has  single-phase  turbulent  forced 
convection  heat  transfer.   The  correlations  associated  with 
this  phenomenon  are  accurate  enough  to  justify  confidence  in 
their  use.   The  fouling  of  the  inside  of  the  tube  is  poorly 
understood,  and  varies  greatly,  so  that  a  necessary  design 
allowance  is  made  and  a  performance  margin  created.   The 
heat  transfer  through  the  tube  wall  itself  is  by  conduction, 
a  classical  problem  in  cylindrical  coordinates,  and  can  be 
calculated  as  accurately  as  the  properties  of  the  tube  material 
and  the  surrounding  heat  transfer  coefficients  are  known. 

Aside  from  a  greater  knowledge  of  coolant-side  fouling 
and  how  to  control  it,  it  is  largely  on  the  steam  side  of  the 
tube  that  advances  of  knowledge  can  be  reasonably  expected. 
Here  there  are  three  important  phenomena.   The  first  is  the 
transfer  of  heat  through  a  thin   layer  of  condensate  from 
a  vapor,  either  quiescent  or  flowing,  which  causes  that  vapor 
to  lose  its  latent  heat  of  vaporization  and  join  the  film  of 


condensate  on  the  tube.   As  condensing  particles  join  this 
liquid  layer  they  retain  some,  none,  or  all  of  the  momentum 
they  possessed  as  a  vapor.   The  dynamic  interaction  of  the 
vapor  and  the  liquid  under  both  forced  flow  and  quiescent 
vapor  conditions  presents  a  formidable  problem  to  the 
researcher.   A  second  phenomenon  on  the  steam  side  of  the 
tube  is  the  presence  of  noncondensi bl e  gases  which  complicates 
and  usually  degrades  the  performance  of  the  condenser.   The 
third  phenomenon  is  the  inundation  of  a  condenser  tube  by 
"condensate  rain,"  or  condensate  condensed  elsewhere  in  the 
condenser  which,  as  it  makes  its  way  to  the  hotwell,  impinges 
upon  the  condenser  tube. 

C.   IMPROVING  THE  HEAT  TRANSFER 

If  the  rate  of  heat  transfer  for  a  given  area  of  condensing 
surface  were  to  be  increased,  then  the  total  size  of  a  con- 
denser could  be  reduced.  The  result  would  be  savings  in  space, 
weight,  and  possibly  cost.   Various  methods  of  increasing 
the  heat  transfer  in  a  naval  condenser  have  been  proposed. 
Among  these  are: 

1 .   Reducing  the  Internal  Fouling 

The  reduction  of  internal  fouling  would  certainly 
increase  the  rate  of  heat  transfer.   While  various  tube  metals 
and  chemical  treatments  to  the  coolant  side  of  the  condenser 
can  reduce  fouling,  operational  and  other  considerations 
have  prevented  their  use  in  naval  condensers.   Aside  from 
mechanical  cleaning,  no  acceptable  method  of  significantly 


reducing  marine  fouling  in  a  naval  condenser  has  presented 
i  tsel f . 

2 .  Eliminating  Noncondens i bl e  Gases 

The  elimination  or  reduction  of  noncondensi ble  gases 
in  the  condensing  steam  would  be  expected  to  increase  the 
heat  transfer  rate.   It  is  a  fact  of  realistic  operation  in 
a  marine  environment  that  noncondensi bl e  gases  cannot  be 
eliminated.   There  is  room  for  improvement  in  their  continual 
removal  from  the  condenser  by  air  ejector  devices.   As  long 
as  complicated  systems  with  mechanical  seals  are  used  in  the 
shipboard  environment,  some  amount  of  noncondensi bl e  gases 
will  be  in  the  condenser. 

3 .  Enhancing  Heat  Transfer 

a.  Shaping  of  Tubes 

Various  ingenious  shapes  of  tubes  have  been  de- 
monstrated to  increase  heat  transfer  rate  for  a  given  amount 
of  tube  surface  area.   Complications  arise  from  these  tubes, 
however,  not  the  least  being  a  sharp  increase  in  the  pumping 
power  required  to  maintain  an  adequate  coolant  flow  rate. 

b.  Promoting  Dropwise  Condensation 

The  use  of  various  promoters  upon  the  tube  outer 
surface  to  prevent  the  forming  of  a  liquid  film,  thus  keeping 
the  condensate  in  droplet  form,  has  been  proposed  and  demon- 
strated by  many  authors.   Since  it  is  not  the  intention  of 
this  thesis  to  examine  the  effect  of  this  phenomenon,  suffice 
it  to  say  that  until  now  there  remains  a  question  as  to  whether 
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any  promoter  has  demonstrated  the  robustness  to  withstand 
continued,  long-term  use. 

4.  Eliminating  or  Reducing  Inundation 

Since  condensers  are  formed  with  banks  of  tubes,  it 
is  inevitable  that  condensate  draining  from  one  tube  should 
impinge  upon  some  other  lower  tube.   The  effect  of  this  in 
filmwise  condensation  is  normally  to  thicken  the  condensate 
film  and  to  reduce  the  rate  of  heat  transfer.   The  proper 
placement  of  collection  trays  or  baffles  is  of  great  impor- 
tance in  channelling  falling  condensate  away  from  lower  tubes. 
The  number  and  placement  of  these  collection  trays  is  not 
examined  in  this  thesis. 

5 .  Utilizing  High  Vapor  Velocities 

In  many  regions  of  a  naval  condenser,  vapor  velocities 
of  300  ft/sec  or  more  may  occur.   These  high  velocities  create 
a  forced  flow  situation  on  the  steam  side,  thin  the  film  of 
condensate,  and  generally  increase  the  heat  transfer  coeffi- 
cients on  the  steam  side.   The  higher  velocities  also  cause 
complications,  one  of  which  is  the  separation  of  the  vapor 
boundary  layer  causing  the  condensate  film  to  thicken  suddenly 
The  heat  transfer  in  the  region  of  the  tube  which  is  after 
the  vapor  phase  separation  point  is  relatively  unknown. 


D.   REVIEW  OF  GENERAL  RESEARCH  IN  CONDENSATION 

In  analyzing  the  condensation  of  steam  on  tubes,  the 
cornerstone  of  research  was  laid  by  Nusselt  [4]  who  analyzed 


a  quiescent  vapor  condensing  on  a  vertical  flat  plate  and  on 
a  horizontal  cylinder.   Nusselt's  analysis  allowed  for  no 
interfacial  shear  stress  between  the  flowing  condensate  and 
the  condensing  vapor.   He  also  used  a  constant  tube  wall  sur- 
face temperature . 

Rohsenow,  et  al . ,  [5]  later  expanded  this  work  to  include 
the  effect  of  vapor  shear  at  the  vapor-liquid  interface  upon 
the  condensation  rates.   This  analysis  was  for  a  vertical 
flat  plate  of  uniform  temperature  with  either  turbulent  or 
laminar  flow  of  condensate. 

Sparrow  and  Gregg  [6]  conducted  the  first  analysis  of  sur- 
face condensation  using  boundary  layer  theory.   They  included 
momentum  terms  in  the  condensate  film  analysis  to  account  for 
acceleration  of  the  flowing  condensate.   They  showed  that  the 
momentum  terms  were  insignificant  under  most  conditions. 

Koh,  et  al .  ,  [7]  did  further  work  allowing  the  falling 
condensate  to  drag  the  vapor  along  with  it,  while  still  using 
a  constant  wall  surface  temperature. 

Shekriladze  and  Gomelauri  [8]  expanded  upon  the  earlier 
work  of  Rohsenow,  et  al . ,  [5],  analyzing  the  effect  of  a 
flowing  vapor  upon  laminar  filmwise  condensation  in  the 
absence  of  significant  gravity  effects.   Their  model  used  a 
constant  wall  surface  temperature. 

Denney  and  Mills  [9]  expanded  this  area  of  analysis  to 
include  the  influence  of  gravity  upon  the  condensate  film, 
holding  to  a  constant  wall  surface  temperature. 
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E.   RECENT  WORK  CONCERNING  CONDENSATION  ON  TUBES 
1  .   Vapor-Condensate  Interfacial  Shear  Models 

One  of  the  areas  of  continual  interest  in  condensation 
of  steam  on  a  tube  is  the  question  of  the  effect  of  the  shear 
or  drag  of  the  condensing  vapor  upon  the  flowing  condensate. 
Theories  which  attempt  to  quantify  the  interaction  abound,  but 
recent  ones  take  one  of  two  general  forms. 

The  first  type  of  model  says  that  total  shear  is  the 
result  of  the  summation  of  dry  friction  shear  of  the  vapor 
upon  a  bare  tube  and  a  shear  caused  by  the  condensing  vapor 
retaining  some,  none,  or  all  of  its  vapor  phase  momentum.   The 
deceleration  of  this  vapor  as  it  joins  the  liquid  layer  causes 
this  shear.   Various  authors  [8,  10,  11,  12,  13,]  subscribe 
to  some  combination  of  these  two  shears. 

Other  researchers  [14,  15]  subscribe  to  a  theory  which 
reaches  more  deeply  into  the  molecular  level  interactions  and 
in  reality  appears  to  be  a  more  minute  examination  of  the 
same  idea.   Using  a  Reynolds  Flux  concept,  these  authors  state 
that  for  a  flowing  fluid  there  is  a  transport  of  momentum  in 
a  direction  normal  to  the  mass  flow  which  is  on  the  molecular 
level  in  laminar  flow  but  involves  larger  amounts  of  fluid  in 
turbulent  flow.   This  transport  is  a  function  of  temperature 
in  laminar  flow  and  of  temperature  and  free  stream  turbulence 
in  turbulent  flow.   As  one  of  these  cross-flowing  particles 
approaches  a  solid  wall,  it  is  decelerated  to  a  standstill. 
As  it  rebounds  from  the  wall,  it  is  reaccel erated .   Thus, 
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an  overall  rate  of  condensation  can  be  represented  by  super- 
inposing  a  velocity  perpendicular  to  the  wall,  similar  to 
suction  through  a  dry  wall. 

The  latter  theory  includes  within  it  a  relationship 
between  the  rate  of  condensation  and  temperature  drop  that 
must  necessarily  occur  in  the  process  [14].   The  former 
frictional  shear  force  theories  commonly  assume  that  the 
outer  surface  of  the  condensate  film  is  at  the  steam 
saturation  temperature. 

2 .   Wall  Temperature  Models 

Fujii,  et  al.,  [16]  display  from  their  experimental 
data  that  tube  temperature  varies  around  the  outside  surface 
of  the  tube.   In  their  calculations,  however,  they  use  averaged 
values  for  tube  wall  temperature  and  for  the  temperature-depen- 
dent heat  transfer  coefficients  to  calculate  heat  transfer. 

Nobbs  and  Mayhew  [17]  show  that  for  steam  flowing 

vertically  downward,  a  variable  temperature  analysis  results 
in  a  marked  overall  change  in  heat  transfer.   The  state  that 
earlier  isothermal  tube  wall  surface  models  were  optimistic 
in  their  predicted  heat  transfer. 


Since  the  writing  of  this  theis,  new  results  have  been 

obtained  by  Fujii  and  his  coworkers.  The  reader  is  referred 

to  the  Proceedings  of  the  Workshop  on  Modern  Developments  in 

Mari  ne  Condensers ,  Naval  Postgraduate  School,  26-28  March,  1980 
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Nicol  and  Wallace  [13]  show  that  local  heat  transfer 
coefficients  are  indeed  a  function  of  local  tube  surface 
temperature  and  recognize  the  importance  of  the  circumferential 
path  in  heat  conduction.   In   one  of  the  first  analyses  to 
include  this  path  of  heat  transfer,  these  authors  apply  a 
relaxation  technique  and  iterative  scheme,  changing  the 
temperature  distribution  and  then  the  local  heat  transfer 
coefficients.   It  is  interesting  to  note  that  their  results 
show  that  when  the  circumferential  heat  conduction  path  is 
included,  the  average  heat  transfer  coefficient  rises  about 
twenty -five  percent. 

Nicol,  et  al  .  ,  [18]  analyze  steam  crossflow  upon  condenser 
tube  (Fig.  1)  using  a  formulation  wherein  local  tube  surface 
temperatures  are  used  to  calculate  local  heat  transfer  co- 
efficients.  Thier  results  show  that  the  variable  wall 
temperature  formulation  predicts  less  heat  transfer  then  older 
isothermal  wall  models.   The  effect  of  circumferential  con- 
duction around  the  tube  wall  is  not  included  in  their  analysis. 
A  close  look  at  Fig.  7  of  their  work,  partially  reproduced  here 
as  Fig.  2,  shows  a  theoretical  wall  temperature  profile  and 
a  few  measurements  of  surface  temperature  around  the  tube  wall. 
In  their  discussion  of  this  particular  figure  they  state: 

"The  profiles  are  seen  to  agree  quite  well,  the  main 
diferrences  being  in  the  rounding  off  of  the  more  extreme 
peaks  of  the  theoretical  profile  because  of  the  circumferential 
conduction  in  the  tube  wall,  which  is  of  course  more  pronounced 
in  those  regions  where  the  temperature  gradient  is  high." 

T3 


Nicol,  et  al.,  [18]  also  present  theoretical  curves  of 
local  heat  transfer  coefficients,  local  heat  flux,  and  local 
tube  surface  temperature.   These  cases  are  presented  for  both 
an  older  isothermal  wall  model  and  an  anisothermal  wall  model 
where  heat  can  flow  only  radially  in  the  tube  wall. 

Fujii,  et  al.,  [19]  use  the  results  of  Nicol,  et  al., 
[18]  for  variable  wall  surface  temperature,  coupled  with 
earlier  work  [20]  and  develop  another  model  based  upon  the 
assumption  of  uniform  heat  flux  around  the  tube.   These 
authors  also  recognize  the  path  of  circumferential  heat  con- 
duction as  an  important  area  for  further  work. 

F.    Purpose  of  Research 

It  is  in  this  final  area  that  the  author  has  analyzed, 
with  the  computer,  a  single  tube  of  a  condenser  and  has 
included  in  his  analysis  the  circumferential  conduction  of 
heat  in  the  tube  wall,  a  phenomenon  not  customarily  included 
in  condenser  heat  transfer  analysis. 
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II.   A  SPECIFIC  PROBLEM  FOR  STUDY 

A.  THE  SINGLE  TUBE  PROBLEM 

A  review  of  prior  research  has  indicated  ideas  and  trends. 
One  basic  idea  was  that  local  heat  transfer  coefficients  de- 
pend upon  local  temperature  differences  and  that  the  shearing 
force  depends  upon  these  heat  transfer  coefficients.   The 
trend  was  away  from  older  isothermal  wall  models  to  newer 
variable  wall  temperature  (ani sothermal )  models.   The  logical 
continuation  of  the  trend  would  indicate  that  a  wall  fully 
capable  of  two  dimensional  heat  conduction  would  have  been 
better  still. 

Since  actual  tube  walls  can  conduct  heat  in  the  circum- 
ferential direction,  then  a  model  which  could  account  for 
this  would  be  an  improvement  over  previous  models.   In  order 
to  test  the  plethora  of  shear  force  and  condensation  rate 
theories,  a  model  including  the  effect  of  circumferential 
conduction  is  desirable,  and  may  prove  to  be  essential. 

B.  THE  PROBLEM  STATEMENT 

The    problem  was    posed    for   the    upper   half   of   a    cylindrical 
condenser    tube   with    a    cross-flowing    vapor,    that    is,    where    the 
direction    of    the    gravitational    force    and    the    direction    of 
vapor    flow   where    perpendicular,    (See    Fig.    la).      Symmetry   with 
the    other   half   of   the    tube   was    assumed,    allowing    the    simplifi- 
cation   that    the    ends    of   the    half    tube,    where    it   would    normally 
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join    the    other    half,    were    adiabatic   walls.      The    problem    for 
solution   was    set    up    so    that   9,    the    local    temperature    differ- 
ence,   was    defined    thusly: 

e(r.f)    =   Ts    -    T(r,o)  (1) 

While    T      may    vary   with    pressure    changes    around    the    tube,    it 
is    assumed    constant   herein.      The    governing   equation    for   heat 
transfer   within    the    tube   wall    is    the    Laplace    Equation: 

1       39     .    1       329 
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with    the    following   boundary    conditions: 
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On    the    inside    of   the    tube    turbulent    forced    convection   was 
assumed    and    the    Di ttus-Boel ter    correlation    [21]   was    used    to 
provide    a   heat   transfer   coefficient. 

On    the   outside    of    the    tube,    the    condensate    flow   was    taken 
to    be    laminar,    and    heat    transfer    through    the    thin    condensate 
layer   was    assumed    to    be    purely    by    conduction.      For   such    con- 
ditions   the    film   heat    transfer    coefficient   is    given    by: 


(3) 


The  film  thickness  5  was  calculated  from  the  formulation  due 
to  Nicol,  et  al.,  [18]  (see  also  Appendix  A).  This  analysis 
yields    the    following    first-order   differential    equation: 


3*rok*9o 
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The    shear    term,    x    ,    used    by    Nicol,    et    al.,    [18]    was    one 
given    by   Schlicting    [22]    and    is    of    the    following    form: 

t      =    ,,  ,(6.973$    -    2.732<J,3    +    0.292(f)5    -    0.0183<{.7 


fuvropv^ 

V     uv      I      +   0. 


000043cj>9    -    0.00011  54.11)  .  (5) 

For    the    sake   of   comparison    this    author    used    the    same    shear 
stress    model.      This    model    calculates    vapor    boundary    layer 
separation    at    about    108    degrees    from    the    forward    stagnation 
point.      Thereafter,    the    shear    stress    was    set    to    zero. 

In   order    to    integrate    Eq.    (4)    a    starting    value   of   con- 
densate   film   thickness    at    the    stagnation    point   was    needed. 
Because    high    velocity    steam    flowing    horizontally    normal    to    a 
horizontal    tube   was    to    be    studied,   -7—  was    set    equal    to    zero 
for    the    stagnation    point    (<j>    =0),    and    an    approximate    starting 
value    for    the   condensate    film    thickness,    <S ,    was    found    from 
Eqs  .    (4)    and    (5) . 

Once    this    starting    value   was    in    hand,    Eq.    (4)    was    inte- 
grated,   in    its    entirety,    around    the    tube    using    an    estimate   of 
the    tube    surface    temperature    difference,    [9      =   e(r    ,    *)].      The 
integration    resulted    in    the    local    film   thickness    around    the 
tube   wall.      Local    heat   transfer    coefficients    for    the    steam 
side   were    then   calculated    according    to    Eq .    (3). 

Using    these    local    heat    transfer   coefficients    Eq.    (2)    could 
then    be    integrated,    resulting    in    a    temperature    distribution 
throughout    the    solid    tube   wall.      With    a    new    set   of   tube-wall 
outer   surface    temperatures    in    hand,    a    new   set   of   local    heat 
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transfer  coefficients  could  be  calculated  and  applied  for 
another  iteration  of  the  Laplace  equation.   If  temperatures 
converged  and  heat  transfer  balanced,  then  a  solution  was  in 
hand  . 

C.    THE  GALERKIN  FINITE  ELEMENT  METHOD 

The  purpose  of  this  application  of  the  Galerkin  Finite 
Element  method  is  to  solve  a  steady  state  two  dimensioanl 
heat  conduction  problem  by  converting  the  governing  partial 
differential  equation,  the  Laplace  equations  in  this  case, 
to  a  series  of  algebraic  equations.   Each  of  these  algebraic 
equations  represents  the  partial  differential  equation  over 
a  samll  sub-region  of  the  total  domain  called  an  element. 
The  result  of  the  simultaneous  solution  of  these  algebraic 
equations  is  to  drive  the  integral  square  error  of  the 
approximate  solution  to  zero  over  each  element.   Elements  may 
be  of  different  sizes  and  may  be  clustered  in  areas  where 
greater  accuracy  is  desired. 

The  Galerkin  Finite  Element  method  differs  from  the 
Galerkin  method  in  that,  since  the  algebraic  equation  repre- 
senting the  approximate  solution  exists  over  only  one  small 
element,  the  resultant  matrix  of  algebraic  equations  for 
solution  takes  on  a  banded  form.   That  is,  it  has  zero  upper 
and  lower  triangles  of  significant  size.   The  resultant 
savings  in  computer  storage  space  and  computational  time  by 
using  special  matrix  equation  solvers  is  significant. 
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In  this  discreet  approximation  method  the  boundary  con- 
ditions are  applied  locally  to  each  element  on  the  boundary 
of  the  domain.   In  this  case  linear  triangular  elements 
were  used.   The  algebraic  approximations  to  the  solution, 
called  basis  or  shape  functions,  were  linear.   The  resultant 
approximate  solution  amounts  to  approximating  a  smooth  three 
dimensional  surface  by  a  mosaic  containing  a  large  number  of 
small  triangles.   Triangles  size  varied,  small  ones  being 
used  in  such  areas  of  interest  as  the  forward  stagnation 
point  and  anywhere  a  steep  temperature  gradient  could  be 
expected  to  exist.   Larger  elements  were  used  in  areas  where 
the  temperature  profile  was  considered  to  be  relatively 
unchanging.   As  solutions  progressed,  refinement  of  the 
triangular  meshes  highlighted  and  concentrated  upon  areas  of 
interest. 

For  greater  detail  on  both  the  Galerkin  method  and  the 
Galerkin  finite  element  method,  the  interested  reader  is 
referred  to  Ozisik  [23]  and  Fairweather  [24]. 

The  formulation  and  computer  program  were  tested  using 
classical  two-dimensional  problems  of  heat  conduction  from 
Carslaw  and  Jaeger  [25].   The  first  of  these  was  a  rectangle 
with  one  side  at  a  given  temperature  and  the  other  three  at 
a  specified  lower  temperature.   The  second  was  a  rectangle 
with  two  adjoing  sides  adiabatic,  a  third  side  at  a  given 
temperature,  and  the  fourth  having  a  given  heat  transfer 


19 


coefficient  to  a  medium  at  some  lower  temperature.  In  both 
cases  the  finite  element  solutions  converged  to  exact  solutions 

A  further  test  of  the  program  and  formulation  was  the 
problem  of  a  rectangle  with  two  opposite  sides  adiabatic,  a 
theird  side  at  a  given  temperature,  and  the  fourth  side  having 
a  local  heat  transfer  coefficient  to  some  medium  at  a  lower 
temperature.   This  local  heat  transfer  coefficient  was 
dependent  upon  both  the  local  surface  temperature  and  the 
position  along  the  wall.   Since  this  problem  cannot  be  compared 
to  a  classical  solution,  any  solution  could  only  be  given 
credence  by  detailed  examination.   It  was  seen  that  the 
solution  converged  for  smaller  and  smaller  triangle  sizes, 
that  plots  of  isotherms  internal  to  the  rectangle  were  smooth 
and  not  obviously  in  error,  and  that  heat  transfer  into  and 
out  of  the  rectangle  was  in  balance.   While  certainly  no 
proof,  these  comparisons  led  to  confidence  in  the  program. 

D.    APPLICATION  OF  THE  INNER  BOUNDARY  CONDITION 

For  the  inside  of  the  tube,  the  boundary  condition  of  con- 
vection heat  transfer,  Eq.  (2),  can  be  simple  represented  by 
noting  that  the  inside  film  coefficient,  h.,  is  given  by  the 
Di ttus-Boel ter  correlation  [21],  and  is  a  function  of  coolant 
bulk  temperature  and  flow  rate.   Thus,  since  the  thermal 
resistance  at  his  interface  is  independent  of  local  tube  wall 
temperature,  the  resistance  may  be  represented  by  a  fictitious 
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thickness    of    insulating   material    according    to    the    relationship 


R    = 


1 


hiA 


ri    ln    (r/1    tf} 


KfA 


(6) 


From    this,    the    thickness    of   the    fictitious    ring    tf,    was    deter- 
mined   and    the    equivalent   boundary    condition,    that    the    coolant 
temperature    existed    at    the    inside    of    the    fictitious    ring,    was 
i  mposed . 

E.   SOLUTION  PROCEDURE 

On  the  outside  of  the  tube,  the  local  heat  transfer  co- 
efficients were  computed  all  around  the  tube  for  the  current 
value  of  the  calculated  local  outside  wall  temperature.   These 
local  heat  transfer  coefficients  were  applied  to  elements 
having  sides  on  the  outer  tube  wall. 

In  order  to  converge  to  the  proper  wall  temperatures  and 
temperature-dependent  heat  transfer  coefficients,  an  iterative 
scheme  was  followed  wherein  a  new  set  of  local  heat  transfer 
coefficients  for  the  outside  of  the  tube  were  computed  from 
new  solution  temperatures.   The  process  was  repeated  until 
temperatures  at  all  nodal  points  changed  less  than  0.02  degrees 
Celsius  for  consecutive  iterations. 

Having  stabilized  the  local  temperatures  throughout  the 
tube  wall,  a  check  was  made  to  insure  that  the  heat  transfer 
balanced.   The  total  heat  deposited  upon  the  condenser  tube 
outer  wall  by  the  condensing  steam  was  compared  to  the  total 
heat  conducted  into  the  tube  wall  at  the  outer  surface  and  to 
the  total  heat  conducted  from  the  tube  at  its  inner  surface. 
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F.   COMPARISON  WITH  PREVIOUS  SOLUTIONS 

The    aforementioned    paper   of   Nicol,    et    al . ,    [18]    provided 
fertile    ground    for    comparison,    but   only    if   the    test   conditions 
could    be    simulated.      These    authors    were    able    to    achieve    a 
given    average    tube    outer   surface    temperature    by    varying    the 
flowrate    of    the    coolant    and    hence    the    internal    heat   transfer 
coefficient.      In    the    finite    element    formulation    this    condi- 
tion   is    obtained    by    varying    the    thickness    of   the    aforementioned 
fictitious    ring.      In    order    to    preclude    significant   circum- 
ferential   conduction    in    the    fictitious    ring,    the    circumferen- 
tial   thermal    conductivity    of   the    ficticious    ring   was    chosen 
to    be    one    one-hundredth    of   the    thermal    conductivity   of   the 
tube   wal 1 . 
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Ill .   RESULTS 

A.   TUBE  WALL  TEMPERATURE 

The  tube  wall  temperature  predictions  of  the  author's 
finite  element  model  are  shown  in  Fig.  (3),  and  are  seen  to 
support  the  prediction  of  Nicol  et  al . ,  [18]  that  circum- 
ferential heat  conduction  would  substantially  round  off  the 
theoretical  temperature  profile  obtained  by  considering 
radial  heat  conduction  only.   It  is  also  worthwhile  to  note 
that  the  maximum  and  minimum  temperatures  are  less  extreme 
than  those  predicted  by  the  radial  heat  flow  model  of  Nicol 
et  al.,  [18]  and  generally  agree  more  closely  with  their 
experimental  results.   While  the  actual  differences  in 
temperature  appear  small,  near  the  separation  point  of  the 
vapor  flow  these  differences  approach  thirty  percent. 

Nicol  et  al . ,  [18]  also  show  a  theoretical  tube  wall 
surface  temperature  profile  for  another  average  wall  tempera- 
ture.  This  profile  is  reproduced  in  Fig.  4  along  with  a 
theoretical  profile  developed  from  the  author's  model.   As 
before,  the  steep  gradients  are  rounded  off  and  the  overall 
temperature  profile  is  less  extreme. 

A  plot  of  the  isotherms  within  the  tube  wall  is  quite 
revealing  and  is  shown  in  Fig.  (5).   The  direction  of  vapor 
flow  and  point  of  vapor  boundary  layer  separation,  as  indicated 
by  the  thickening  of  the  condensate  layer,  are  shown  on  the 
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figure.   The  isotherms  are  seen  to  start  close  together  near 
the  vapor  stagnation  point  and  to  spread,  reflecting  lower 
rates  of  heat  transfer,  progressively  around  the  tube.   In 
the  vicinity  of  separation  of  the  vapor  boundary  layer,  the 
isotherms  are  seen  to  deviate  significantly  from  their  con- 
centric upstream  orientation.   This  reflects  a  significant 
amount  of  heat  transfer  ci rcumferenti al ly  around  the  tube  in 
this  general  vicinity.   The  much  wider  spacing  of  the  isotherms 
on  the  back  part  of  the  tube  indicates  that  according  to  this 
model  it  is  an  area  of  relatively  low  heat  transfer  rate. 
1  .   Local  Heat  Transfer  Coefficients 

Figure  6  shows  a  comparison  of  local  heat  transfer 
coefficients  around  the  tube  wall.   The  isothermal  wall  and 
radial  heat  flow  local  transfer  coefficients  of  Nicol  et  al  .  , 
are  partially  reproduced.   The  results  of  the  finite  element 
model  are  seen  to  lie  generally  between  the  other  two  models 
with  the  exception  of  the  vapor-separated  region.   In  the 
region  prior  to  separation,  as  seen  in  Fig.  4,  the  wall 
temperature  difference  of  the  finite  element  model  lies  be- 
tween the  isothermal  wall  case  and  radial  heat  flow  case. 
From  Eq.  (4)  at  the  stagnation  point,  the  initial  value  of 
the  film  thickness  is  proportional  to  9   /3.   Therefore,  a 
larger  temperature  difference  will  create  a  larger  initial 
film  thickness,  and  therefore  a  smaller  heat  transfer 
coefficient,  as  given  by  Eq.  (3).   It  is  not  unexpected  that 
this  formulation  and  solution  should  predict  a  greater  heat 
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transfer    coefficient    in    the    separated    region    than    the    iso- 
thermal   wall    model.      Since    there   was    a    higher   wall    temperature 
and    therefore    less    steam   condensed    in    the    unseparated   portion 
of    the    tube    in    the    present   model,    there    is    a    thinner    conden- 
sate   layer    in    the    after   section    of   the    tube.       In    the    separated 
region    the    greater   heat    transfer    coefficient    predicted    by 
the    author's    model    over    the    radial    heat    flow   model    is    not 
fully   explained,    but    is    no-doubt   due,    in    part,    to    the    assumption 
of   symmetry    in    the    present    formulation. 
2 .      Local    Heat    Flux 

In    still    a    third   examination    under    these    same    condi- 
tions,   Nicol    et   al .    in    their    Fig.    6    present    local    heat    flux 
for   both    the    isothermal    wall    and    radial    heat    flow   models. 
These    profiles    are    partially    reproduced    here    in    Fig.    7    along 
with    the    author's    results   which    again,    as    expected,    lie    between 
the    two    previous    models    except    in    the    separated-vapor    region 
of   the    tube.       It    is    interesting    to    note    that    the    inclusion 
of   circumferential    effects    generates    a    heat    flux   which    is 
much    closer    to    being    uniform   than    that   which    is    obtained   with 
the    isothermal    wall    model. 
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IV.   CONCLUSIONS 

The  circumferential  path  of  heat  conduction  has  been 
successfully  included  in  the  solution  of  laminar  filmwise 
condensation  of  high  velocity  steam  on  condenser  tubes.   The 
author's  anisothermal  fully-conducting  wall  model  yields 
closer  agreement  to  previously  measured  tube  wall  surface 
temperatures  than  a  previous  anisothermal  wall  model  which 
allowed  only  radial  heat  flow.   The  maximum  and  minimum  values 
of  wall  temperature  were  less  extreme  for  the  anisothermal 
fully  conducting  wall  than  for  the  anisothermal  wall  with 
only  radial  heat  flow. 

The  addition  of  a  circumferential  path  for  heat  flow 
results  in  a  complicated  interplay  between  this  heat  flow 
and  heat  being  deposited  locally  by  condensation  and  leaving 
in  the  internal  cooling  water.   The  addition  of  vapor  velocity 
effects  on  the  outside  surface  of  the  film,  and  especially 
the  sudden  changes  arising  from  separation  of  the  vapor  flow 
result  in  a  highly  complex  sequence  of  thermal  and  fluid- 
dynamical  interactions.   The  methods  developed  herein, 
though  s ti 1 1  sub ject  to  comprehensive  experimental  verification, 
are  adequate  to  predict  these  complex  phenomena. 
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V.   RECOMMENDATIONS 

Now   that   a   model    for    two-dimensional    heat    conduction   which 
includes    temperature-dependent    local    heat    transfer    coefficients 
has    been    developed    and    demonstrated,    the    use    of    this    model    for 
further   study    is    suggested.      Areas    thought    fruitful    for    further 
examination    are    both    the    shear   stress    models    and    the    tempera- 
ture-dependent  condensation    rate    models.      For   shear    stress 
models    the    separated    vapor    region   or    trailing    area   must   be 
given    closer   examination    since    it    accounts    for   about    forty 
percent   of    the    surface    area    of    the    tube    and    the    fluid    flow 
characteristics    in    this    region    remain    relatively    unknown.      In 
the    study    of    temperature-dependent    condensation    rate    models, 
further   examination    of    the    effect   of    the    velocity    of    the 
flowing    vapor   and    local    temperature    is    expected    to    be    fruitful. 

The    vigorous    exercise    of    this    model    through    a    broad    range 
of   steam-side    and    coolant   side    conditions    is    suggested. 

The    model    of   condensation   with    constant   heat    flux    pro- 
posed   by    several    authors    [18,    20,    26],    should    be    re-examined 
since    circumferential    heat   conduction    in    the    tube   wall,    when 
included    in    theanalysis,   seems    to   yield    a    nearly    constant 
heat    flux   distribution. 
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APPENDIX  A 

Local  Rate  of  Change  of  Condensate  Film  Thickness 

A  force  balance  on  the  shaded  element  of  Fig.  lb  yields 
the  following: 

Tvrod(f)  "  Tirod<^    '    S^-Py)005*  r0  (5-y)  d*  =  ° 


so  that 


3U 


To  =  ^T7  =  Tu  -  (p„-p„)9  cos*  (5-y) 
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Integrati  ng 
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The    average    velocity,    u    is: 
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.    or 
(8) 


The  mass  flow  rate  of  condensate  around  the  tube  is: 

m„  =  p„$u 
and  ,  wi  th  Eq .  ( 8) , 

m„  =  -a 5Z  -  *■* *-  cos*  5° 

The    rate    of   change    of   mass    flow    rate   with    respect    to    angle    is 


<%     =  P^2  <*%   +  V/     d£  +  gp3.(pjTpv)63   sin* 
d*  2y       d*  u         d*  3y 


_   gp£(pH-pv)    cos*   S2   dfi_ 
u  d* 
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Since    drfl/d<|>    can    also    be    equated    to    the    amount    of    steam 
condensing,    the    following    relationship    applies: 


dm 
d$~ 


r    k„9  - 

0     i     0 


5h' 


fg 


where  h1  .   is  the  latent  heat  of  vaporization  corrected  for 

fg 

subcooling.      By   making    this    substitution    and    solving    for 
d6/dct>,    the    relationship    given    in    Eq.    (4)    is    obtained. 
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Figure  1 .  Schematic  diagrams  of  physical  model 
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Figure  2.     Predicted  and  measured  T     versus  angle  from  stagnation  point  [18] 
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Figure  3.     Predicted  and  measured  T     versus  angle  from  forward  stagnation 
point. 
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Figure  4.  Predicted  T  versus  angle  from  forward  stagnation  point 
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Figure  6.  Predicted  local  heat  transfer  coefficient  versus  angle  from 
forward  stagnation  point. 
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Figure  7.  Predicted  local  heat  flux  versus  angle  from  forward  stagnation 
point. 


51 


LIST  OF  REFERENCES 


2. 


Standards  for  Steam  Surface  Condensers,  6th  ed.,  Heat 
Exchange  Institute,  1970. 

Standards  of  Tubular  Exchanger  Manufacturers  Association, 
4th  ed.,  Tubular  Exchanger  Manufacturers  Association,  Inc., 
1959. 


3.  Department  of  the  Navy,  Design  Data  Sheet  Section  D0S4601-1, 
Steam  Condensers  ,  15  October  1953. 

4.  Collier,  J.  G.  and  Wall  is,  G.  B.,  Two  Phase  Flow  and  Heat 
Transfer ,  lecture  notes  presented  for  a  summer  course, 
July  24  -  August  4,  1967,  Stanford,  California,  Volume  4, 
pp  1210  -  1215,  1967. 

5.  Rohsenow,  W.  M.,  Webber,  J.  H.,  and  Ling,  A.  T.,  "Effect 
of  Vapor  Velocity  on  Laminar  and  Turbulent  -  Film  Conden- 
sation," Trans.  ASME,  V.  78.  pp  1637  -  43,  1956. 

6.  Sparrow,  E.  M.  and  Gregg,  J.  L.,  "A  Boundary  Layer  Treat- 
ment of  Laminar  Film  Condensation,"  Trans  .  ASME ,  V.  81 C , 
pp  13  -  18,  1959. 

7.  Koh,  J.  C.  Y.,  Sparrow,  E.  M.,  and  Hartnett,  J.  P.,  "The 
Two-Phase  Boundary  Layer  in  Laminar  Film  Condensation," 
Int.  J.  Heat  Mass  Transfer,  V.  2,  pp  69  -  82,  1961. 

8.  Shekriladze,  I.  G.  and  Gomelauri,  V.  I.,  "Theoretical 
Study  of  Laminar  Film  Condensation  of  Flowing  Vapor,"  Int . 
J.  of  Heat  Mass  Transfer,  V.  9,  pp  581  -  591,  1966. 


9.   Denny,  V.  E.  and  Mills,  A.  F.,  "Non-Similar  Solutions  for 
L.aminar  Film  Condensation  on  a  Vertical  Surface,"  Int . 
J.  of  Heat  Mass  Transfer,  V.  12,  pp  965  -  979,  1969. 

10.  Mayhew,  Y.  R.,  Griffiths,  D.  J.,  and  Phillips,  J.  W., 
"Effect  of  Vapor  Drag  on  Laminar  Film  Condensation  on  a 
Vertical  Surface,"  Proc.  Instn.  of  Mech.  Engrs.,  V.  180 
(3J) ,  pp  280  -  289,  1965  -  66. 

11.  Mayhew,  Y.  R.  and  Aggarwal  ,  J.  K.,  "Laminar  Film  Conden- 
sation with  Vapor  Drag  on  a  Flat  Surface,"  Int .  J .  of 
Heat  Mass  Transfer,  V.  16,  1944  -  49,  1973. 


52 


12.  Denny,  V.  E.  and  South  III,  V.,  "Effects  of  Forced  Flow, 
Noncondensabl es ,  and  Variable  Properties  on  Film  Conden- 
sation of  Pure  and  Binary  Vapors  at  the  Forward  Stagnation 
Point  of  a  Horizontal  Cylinder,"  Int.  J .  of  Heat  Mass 
Transfer,  V.  15,  pp  2133  -  2142,  1972. 

13.  Nicol,  A.  A.  and  Wallace,  D.  J.,  "Condensation  With 
Appreciable  Vapor  Velocity  and  Variable  Wall  Temperature," 
Steam  Turbine  Condensers,  NEL  Report  No.  619,  pp  27  -  38, 
1976. 

14.  Silver,  R.  S.,  "An  Approach-  to  a  General  Theory  of  Surface 
Condensers,"  Proc.  Instn.  Mech.  Engrs.,  V.  178,  pp  339  - 
376,  1963  -  64. 

15.  Wallis,  G.  B.,  "Use  of  the  Reynolds  Flux  Concept  for 
Analyzing  One-Dimensional  Two  Phase  Flow,  Int .  J .  of  Heat 
Mass  Transfer,  V.  11,  pp  445  -  458  and  459  -  472,  1968. 

16.  Fujii,  T.,  Vehara,  H.  and  Kurata,  C,  "Laminar  Filmwise 
Condensation  of  Flowing  Vapor  on  a  Horizontal  Cylinder," 
Int.  J.  Heat  Mass  Transfer,  V.  15,  235  -  246,  1972. 

17.  Nobbs,  D.  W.  and  Mayhew,  Y.  R.,  "Effect  of  Downward  Vapor 
Velocity  and  Inundation  on  Condensation  Rates  on  Horizon- 
tal Tube  Banks,"  Steam  Turbine  Condensers,  NEL  Report  No. 
619,  pp  39  -  52,  1976. 

18.  Nichol,  A.  A.,  Bryce,  A.,  and  Ahmed,  A.  S.  A.,  Condensation 
of  a  Horizontally  Flowing  Vapor  on  a  Horizontal  Cylinder 
Normal  to  the  Vapor  Stream,  paper  presented  at  Sixth  Inter- 
national Heat  Transfer  Conference,  Toronto,  Canada,  August 

7  -  11  ,  1978. 

19.  Fujii,  T.,  Honda,  H.  and  Oda,  K.,  "Condensation  of  Steam  on 
a  Horizontal  Tube  -  The  Influence  of  Oncoming  Velocity  and 
Thermal  Conduction  at  the  Tube  Wall,". 

20.  Fujii,  T.,  Uehara,  H.,  and  Oda,  K.,  "Filmwise  Condensation 
on  a  Surface  with  Uniform  Heat  Flux  and  Body  Force  Con- 
vection," Heat  Transfer  -  Japanese  Research,  V.  1,  no.  4, 
pp  76  -  83,  1972. 

21.  Holman,  J.  P.,  Heat  Transfer,  p  204,  McGraw  -  Hill,  1976. 

22.  Schlichting,  H.,  Boundary  Layer  Theory,  p  160,  McGraw  - 
Hill  ,  1968. 

23.  Ozisik,  M.  N.,  Boundary  Value  Problems  of  Heat  Conduction, 
pp  338  -  345,  International  Textbook  Co.,  1968. 


53 


24. 

25. 
26. 


Fai  rweather ,  G . . 

Differential  Equations,  pp  21  -  66 

1978. 


Finite  Element  Galerkin  Methods  for 


Marcel  Dekker  ,  I nc , 


Carslaw,  H.  S.  and  Jaeger,  J.  C,  Conduction  of  Heat  in 
Sol  ids ,  pp  166  -  168,  Oxford:  Clarendon  Press,  1959. 

Mills,  A.  F.,  Tan,  C.,  and  Chung,  D.  K.,  "Experimental 
Study  of  Condensation  from  Steam-Air  Mixtures  Flowing 
Over  a  Horizontal  Tube:  Overall  Condensation  Rates," 
Proc.  Fifth  Intl.  Heat  Transfer  Conf.,  pp  20  -  23,  1974 


54 


INITIAL  DISTRIBUTION  LIST 


1  . 


2. 


3. 


5. 


7. 


8. 


Defense    Technical 
Cameron    Station 
Alexandria,    Virginia 


Information   Center 
22314 


No.  Copies 
2 


Library,  Code  0142 

Naval  Postgraduate  School 

Monterey,  California   93940 

Department  Chairman,  Code  69 
Department  of  Mechanical  Engineering 
Naval  Postgraduate  School 
Monterey,  California   93940 

Professor  P.  J.  Marto,  Code  69 
Department  of  Mechanical  Engineering 
Naval  Postgraduate  School 
Monterey,  California   93940 

Assoc.  Professor  R.  H.  Nunn,  Code  69 
Department  of  Mechanical  Engineering 
Naval  Postgraduate  School 
Monterey,  California   93940 

Assoc.  Professor  D.  Salinas,  Code  69 
Department  of  Mechanical  Engineering 
Naval  Postgraduate  School 
Monterey,  California   93940 

Professor  G.  Cantin,  Code  69 
Department  of  Mechanical  Engineering 
Naval  Postgraduate  School 
Monterey,  California   93940 

LT  H.  M.  Holland,  USN 

139  Pelican  Drive 

Groton,  Connecticut   06340 


55 


-  ■- 


V 


1  9§  W 

H6845    Holland 

c.l         The  effect  of  cir- 
cumferential tube  wall 
heat  conduction  upon 
laminar  filmwise  con- L 
densation  on  the  out- 
side of  condenser 
tubes • 

I     is  APPf  «»^  5  7  8  2  P  « 


Thesis  192417 

H6Qk5  Holland 

c.l         The  effect  of  cir- 
cumferential tube  wall 
heat  conduction  upon 
laminar  filmwise  con- 
densation on  the  out- 
side of  condenser 
tubes. 


The  effect  of  circumferential  tube  wa 


3  2768  002  06901  5 

DUDLEY  KNOX  LIBRARY 


I 


